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Abstract 

We develop a straightforward method to compute charge radii and quadrupole moments for 
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J> ' baryons both with and without strangeness, when the number of QCD color charges is N c . The 

00 ' 

minimal assumption of the single-photon exchange ansatz implies that only two operators are 

m ; 

required to describe these baryon observables. Our results are presented so that SU(3) flavor and 
isospin symmetry breaking can be introduced according to any desired specification, although we 
also present results obtained from two patterns suggested by the quark model with gluon exchange 
interactions. The method also permits to extract a number of model-independent relations; a 
r*| , sample is r\jr\ = 3/(N c + 3), independent of SU(3) symmetry breaking. 
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I. INTRODUCTION 



Between the static properties of hadrons, e.g., masses, electromagnetic moments, matter 
and charge radii, and low-energy dynamical properties such as scattering lengths and decay 
rates, a wide variety of accurately measured data lies outside the perturbative range of QCD. 
While techniques such as lattice gauge theory, QCD sum rules, and wide variety of models 
have been developed to study this extremely interesting energy regime, the only known 
perturbative approach for strongly-interacting Yang-Mills theories that holds at all energies 
is the 1/N C expansion [JXJ] , where N c is the number of color charges. 

The baryon sector is particularly amenable to this expansion, since baryons with N c 
colors contain N c valence quarks. Then, each additional quark participating in an interaction 
(specified by an operator with known transformation properties under spin and flavor) brings 
in a factor of the QCD coupling constant a s oc 1/N C due to the requirement of one or more 
gluons to connect this quark to the interaction. An operator in which n quarks interact 
is called an n-body operator; its suppression in powers of 1/N C tends to increase as n 
increases. This "operator method" has recently been used to study a wide variety of baryon 
observables (see Ref. for a recent list of references). The strength of the 1/N C expansion is 
that additional gluons do not spoil this counting, and the only other powers of N c that need 
be taken into account arise from the combinatorics of the quarks. It should also be noted 
that the glue and sea quark pairs in the baryon, not just the valence quarks, are subsumed 
by the operator method ||. 

This work extends the results of two of our recent papers on baryon observables in the 
1/N C expansion, regarding charge radii || and quadrupole moments ||, to baryons with 
nonzero strangeness. These two types of observables are studied here in a single paper 
because the method of calculation is, as seen below, very similar in the two cases. As 
mentioned in the previous works, a treatment of the strange sector was not undertaken at 
the time of their writing because the group theory in the 3-flavor sector is more involved 
and warrants a separate treatment. The operator method calculation specific to the case 
iV c = 3, otherwise known as the general QCD parameterization M was carried out for baryon 
quadrupole moments in ; in that case the calculation can readily be carried out for states 
with or without strangeness since it is relatively simple to perform calculations using the 
full baryon spin-flavor wave functions with only 3 quarks. 
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Here we demonstrate that the solution of the full 3-flavor iV c -quark problem can nonethe- 
less be handled entirely using SU(2) Clebsch-Gordan (CG) coefficients, ultimately because 
the strange states are related to nonstrange states through the SU(2) U- and V^-spin sub- 
groups of SU(3) ||. In the case of the ground-state baryons, i.e., those belonging to the 
large- N c generalization of the SU(6) 56-plet, the total symmetry of the wave function under 
simultaneous exchange of spin and flavor indices simplifies this procedure considerably. 

In this paper we focus only on those aspects of the problem unique to the strange sector; 
the reader is directed to Refs. @ and || for a more thorough discussion of details of the l/N c 
operator method, as well as outlook for the experimental measurement of baryon quadrupole 
moments. The remainder of this paper is organized as follows: In Sec. II we briefly discuss the 
experimental situation regarding measurement of baryon charge radii. Section |TXT| presents 
the details of the 3-flavor calculation; the most important feature of the analysis of || and 
survives, namely, that only one operator appears at the 2-body level and one at the 



3-body level, leading to a large number of constraints between the observables. Section [TV 
presents a sample of the plethora of results obtained from this calculation, e.g., by including 
frequently used patterns of SU(3) flavor symmetry breaking. Section [V] summarizes. The 
Tables and Appendices contain the extensive results of the calculation in forms designed to 
be most useful to interested researchers. 



II. CHARGE RADII: PROSPECTS FOR MEASUREMENT 

The mean-square charge radius r 2 B of a baryon B is defined through the elastic Sachs 
charge form factor F(q 2 ), a function of the photon 4-momentum transfer q 2 often denoted 
by Ge(q 2 ), by the relation 

(2.1) 



B ~ b Q dq 2 
while the total charge Q = F(0) is omitted if Q = 0. 

The operator giving rise to charge radii is a scalar that does not change electric charge 
or strangeness. Thus, it connects only states with the same total J, I3, and Y. However, 
the operator can connect states of different isospin / since the electromagnetic interaction 
does not conserve this quantum number. Each baryon state thus has a charge radius, and 
in addition there exists a r| A transition charge radius. 

At the time of this writing, only the p, n, and S~ charge radii have been measured. 



A charge radius is the first nontrivial moment of a Coulomb monopole (CO) transition 
amplitude, and thus charge radii of hadrons are typically measured through their coupling 
to Coulomb photons in elastic electron scattering. 

In the case of the proton, a dedicated measurement of the elastic electron-proton cross 
section at very low momentum transfer |7J led to the proton root mean-square charge radius 
r p = ^Jr^ = 0.862(12) fm [r 2 = 0.743(21) fm 2 ], a value that is considerably larger than 
the famous dipole result of r p = 0.80 fm obtained by Hofstadter et al. @. A dispersion- 
theoretical analysis of electron scattering data available through 1995 |J gives the result 
r 2 = 0.717(15) fm 2 . However, measurements using the hydrogen Lamb shift as an alternative 



source of information on the proton size tend to produce larger values. For example, Ref. [10 



extracts a value of r p = 0.890(14) fm [r 2 = 0.792(25) fm 2 ], while Ref. [11] extracts a value of 



r 2 = 0.780(25) fm from the experiment of Ref. ||12|| . A reanalysis [13| that includes Coulomb 
and recoil corrections for A > 1 targets produces a result that accommodates both types of 
data: 

r 2 = 0.779(25) fm 2 , (2.2) 
which we use in our analysis. The other available experimental values are 

r 2 = -0.113(3)(4)fm 2 [14], (2.3) 

and the very recently published result 

r|_ = 0.61(12)(9) fm 2 [15]. (2.4) 

The measurement of r|_ in particular suggests the possibility of measuring the charge radii of 
other long-lived strange baryons, the S + , A, H~, S°, and Q~ . In such cases the impediments 
are experimental in nature, particularly the problem of producing a beam of baryons of 
sufficient quality that electron scattering events can be separated from background events. 
But nonetheless one can anticipate such difficulties being overcome in the future. 

For the spin-3/2 resonances and E°, however, the short lifetimes mean that charge radii 
can only be observed in off-shell processes. For example, the A ++ charge radius could be 
extracted in principle through the process ir + p — > n + pe~e + , but in this case one would need 
to use model dependence in separating resonant from continuum n + p scattering, as well as 
isolate the source of the virtual photon 7* — > e + e~ as coming solely from the A. Processes 
with real photons, which for example can be used to measure magnetic moments, cannot 



be utilized to probe the Coulomb transitions from which charge radii are extracted. The 
situation is very similar to that described in Ref. for the baryon quadrupole moments. 
Therefore, our results divide into two categories: Predicitions for observables that may be 
measured in the near-to-medium future, for which the quality of our results can be checked, 
and those that can only be predicted but not measured any time soon. 



III. DETAILS OF THE CALCULATION 



A. Constructing the States 

Naively, one would expect that SU(3) CG coefficients are required to construct a state 
with strangeness. This potential complication can be avoided by noting that the strange 
states are related to the nonstrange states via U- and l/-spin SU(2) subgroups of SU(3). 
In an arbitrary SU(3) representation, one method of calculating matrix elements is to start 
with the nonstrange states and apply the lowering operators U- and V- , imposing along the 
way orthogonality in isospin among otherwise degenerate states such as S° and A, which is 
precisely the same approach as used to obtain flavor eigenstates. 

However, this method does not exploit the full symmetry available to the J p = l/2 + and 
3/2 + ground-state baryon SU(3) multiplets (We do not call them "octet" and "decuplet," 
since the corresponding multiplets for arbitrary iV c are much larger These flavor mul- 

tiplets belong to the large- N c analogue of the SU(6) 56-plet, which is completely symmetric 
under simultaneous exchange of spin and flavor indices. That is to say, the I-, U-, and V- 
spin symmetries are correlated with the SU(2) spin symmetry in a particularly convenient 
way, as we now explore. 

Let N a denote the number of valence quarks of type a in the baryon; then N c = N u + 
Nd + N s and ^3 = (N u — N^)/2. Since the spin-flavor wave function is completely symmetric, 
the spin wave function of quarks of type a within the baryon must be completely symmetric, 
which means that the N a spin-1/2 quarks of type a carry the maximum possible spin value, 
S a = N a /2. This information alone allows to analyze the nonstrange states. In the strange 
baryon case, note that each physical baryon state is still specified by its isospin J, which 
uniquely determines the symmetry of the flavor wave function carried by just u and d quarks. 
Owing to the complete symmetry under spin x flavor, the spin wave function carried by 
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the u and d quarks together must possess exactly the same symmetry properties as the 
corresponding flavor wave function, and hence, the spin S u d = S u + Sd and isospin I have 
the same eigenvalue, S u d = I. 

But now one need only combine the state of ud combined spin S u d = I and isospin 
quantum numbers J, J3 with the symmetrized strange quarks carrying spin S s to obtain the 
complete state with spin eigenvalues J, J3, where J = S u d + S s . To be precise, 

/ 



\jj 3] ii 3 (s u s d s s )) = Yl 



QZ QZ 

\ °ud °s 



E 

Cz qz 



S U Sd 

oz qz 



O, 



ud 



S u ) \Sd Sg) \S S Sg) , 
(3.1) 



where the parentheses denote CG coefficients. Now, in order to compute the matrix elements 
of any particular operator, one need only sandwich it between a bra and ket of the form of 
Eq. (|3.1|) and use the Wigner-Eckart theorem. 



B. The Operator Basis 

In the analysis of both charge radii and quadrupole moments we use the single-photon 
exchange assumption, i.e., that the photon probing these baryon observables couples to 
only one quark line within the baryon. Although physically rather mild, this assumption 
drastically reduces the number of distinct operators that need be considered. (However, 
note that the full basis of operators can certainly be used, as shown for the charge radii || 
or for the magnetic and quadrupole moments |17| .) Using indices i,j, k to indicate quarks 



within the baryon, the most general operator expansion for the charge radius operator out to 
the 3-body level — all that is necessary since 4-, 5-, etc. body operators acting upon physical 
baryons are linearly dependent on those at lower order — reads 

dF(q 2 ) 



-6- 



dq 2 



Nc B Nc C 
= A Y^^i + TF^2 b i°j a i a i + TH E diej f k (Ti-cij, (3.2) 



where the coefficients a, b, c, d, e, / depend only on quark flavor; to obtain the results of 
Ref. ||, set cii = bi = fi = Qi, where Qi is the charge of the ith quark, and c« = di = = 1. 
The corresponding expansion for the quadrupole moment operator reads 

B' C 
Q = TrE & i c j (ZvizVjz - (T ^- (T i) + jH E d 'i e 'j f'k o-i ■ a i , (3-3) 



where the coefficients b' ,d ,d' , e' ', f again depend only on quark flavor, and the results of 
Ref. may be obtained by setting b\ = f- — Qi, d t — d\ — — 1. Note that most models 
use a universal form of SU(3) symmetry breaking for all operators, so that typically the 
corresponding primed and unprimed coefficients are equal. Values used for Qi are discussed 
in Sec. [TV] below. Note that 3-body operators with 3 Pauli matrices are absent [[J due to 
time-reversal symmetry. The use of primed and unprimed coefficients with the same labels 
echoes the similar form of the operators in the two cases; indeed, in a number of models the 



primed and unprimed coefficients are related |18 



One can obtain the most general charge radius and quadrupole moment expressions by 
computing the 12 fundamental matrix elements 

(a/3)® = (S a -Sp), 

(a?)® = (3S^S' p -S a -Sp), (3.4) 

where (a, (3) are the 6 distinct pairs (u,u), (u,d), (u,s), (d,d), (d,s), and (s,s) (Note 
the symmetry under a «-> /3). Then, charge radii or quadrupole moments for baryons in 
the presence of arbitrary SU(3) or isospin symmetry breaking can be computed simply by 
taking an appropriate linear combination of these primitive matrix elements. 

In particular, while the indices i, j, k indicate each of the N c quarks within the baryon, 
let the indices cn,/3, 7 indicate collectively the quarks of a particular flavor u, d, or s. The 
coefficients for the most general operators that appear for the charge radii or quadrupole 
moments can be represented by 3- vectors such as a = (a u , a^, a s ). Thus, matrix elements of 
the 1-body operator for charge radii generalize to: 



a% = N <* a «' ( 3 - 5 ) 

i a 

while those of the 2-body operator may be expanded as 

h Cj ffi-aA = b a c p («/?) (0) - 3 N a b a c a , (3.6) 
and those of the 3-body operator become: 

( E d i e i fktri'Vj) = +4^d„e^ («/?} (0) 

\i^j^k I a,/3 7 

-3 E N * d a e a E Np fp + 6^JV Q ie Q f a , (3.7) 



where N a is the number of quarks of type a (values given in Table |), which satisfies the 
constraints N u + N d + N s = N c and N u - N d = 2I 3 . 

Matrix elements of the 2- and 3-body operators for the quadrupole moments expand as: 

(]T b>c' 3 (3a iz a jz - cT.rcr,)) = 4]T b' a c' p (a^ 2 \ (3.8) 

and 

I /3,7 « 

-4E rf a4(/a + /^)(«/3) (2) - (3.9) 

The two examples of SU(3) symmetry breaking in the quadrupole operators considered in 
the N c = 3 work of Ref. |5j correspond to setting b' a = Q a (m/m a ) n , c' a = e' a = m/m a , 
d' a = (m/m a ) n , and f' a = Q a , where m u = m d = m and m/m s = r represent SU(3) 
symmetry breaking arising from different quark masses. The quadratic (cubic) case in [||] is 
n = 1(2). 

C. Reduction of the Basis 

Even with this decomposition, it is not necessary to compute the entire set of 12 matrix 
elements separately: They are related by a number of simple constraints. First, note that 
Eq. (|3.1|) is sensitive to the exchange of u and d quarks only through the second CG coef- 
ficient, and the factor obtained through this exchange is just (— iy Su+s d- 1 ^ Of course, the 
values of S a , which simply count one-half the number of quarks of flavor a in these baryons, 
remain unchanged from initial to final state. The same is true for I 3 = S u — S d , but the 
total isospin may change to a value I'. One thus finds for an operator O that 

(I'h \0{u <- d)\Ih) = {-I) 1 '' 1 {I' -h P\ I -h). (3.10) 

There are almost enough constraints to allow extraction of all the matrix elements (a(3)^ 
in terms of simple combinations of eigenvalues of the compatible operators of the system. 
To be precise, these are J 2 — > J (J + 1), J z , S 2 a — > (N a /2)(N a /2 + 1) for a = u,d,s, and 
I 2 = S 2 ad = 1(1 + 1). Note that the eigenvalue I3 — S u — S d is not independent. Clearly, 
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(acn)' ' = (N a /2)(N a /2 + 1). The other two constraints are: 

j(j + 1) = (J 2 ) = ((s u + s d + s s ) 2 ) = £(«/?> (0) , 

a,/3 

1(1 + 1) = (S 2 ud ) = ((S u + S d ) 2 ) = (uu)® + 2(ud)M + (drf)(o). (3.11) 

From these constraints it is clear that only one of the two matrix elements (its)' ) or (ds)^ 
need be computed directly from Eq. ( |3.1| ), while the other can be obtained either from the 
inversion rule Eq. (|3.10| ) or from the constraints Eqs. ( |3.11|) . 

As for the matrix elements (a(3)^ 2 \ the most obvious constraint reads: 

(SJ 2 - J 2 ) = (z(S* u + S z d + Sf) 3 - (S u + S d + S s ) 2 ) = EM (2 » (3.12) 

The l.h.s. of this expression is of course simple to compute; in the stretched state (J z = J) in 
which quadrupole moments are computed, it equals J (2 J — 1) for diagonal matrix elements 
and vanishes for transitions. 

Another constraint may be obtained by considering the combination: 

<3&) 2 - S 2 ud ) = (3(J Z - S z s ) 2 - I 2 ) = 3J 2 - 6J Z (S:) + (ss)M + (S 2 S ) - (I 2 ), (3.13) 

whose l.h.s. is just 

(uu)( 2) + 2(ud) {2) + {dd) {2 \ (3.14) 

Note that this constraint requires one to compute also the matrix element (S%), but this 
calculation turns out in purely algebraic terms to be simpler than that of the rank-2 tensors. 
One still requires one more constraint to compute separate values for (us)^ and (ds)^; 
once one of these matrix elements is in hand, the other is obtained by using Eq. ( |3.10| ). The 
required constraint, again using (S^) = (I 2 ), may be obtained from: 

- S U 'J) = 3(S*)J Z - (us)® -\{I 2 + S 2 U - S 2 >, (3.15) 

where the l.h.s. clearly equals: 

+ (us) (3.16) 

In summary, then, the charge radii matrix elements can be computed in complete generality 
using only eigenvalues and the calculation of (its)' ', while the most general quadrupole 



S 



moment matrix elements can be obtained from further explicitly calculating (S*), (S%), 
(uu) {2 \ and (ss)W. 

Another advantage of using the constraints ( |3.11| )-( |3.15| ) is that, by judicious choice of 
which matrix elements to compute explicitly, one may obtain results containing nothing more 
complicated than a 6j symbol. These appear due to the coupling of the quantum numbers 
of expansions for both the bra and ket [Eq. ( |3.1| )] through the operator. For example, in 
the matrix elements of (uu)^ defined by Eq. (|3.4j ), one uses the Wigner-Eckart theorem to 
obtain: 



S U S U 3(5^ 



iz\2 



S' U S':) = 6 SL sM*s^S u (S u + 1)(2S U - 1)(2S U + 3) 



2 S u 



S' 



Cz 

(3.17) 

and the CG coefficient in this expression is linked to the those in the bra and ket of Eq. (p.l|) 



by sums over S* and S' u z . Four appropriately linked CG coefficients produce a 6j symbol, 
and six produce a 9j symbol. In our case, the matrix elements (us)^ 2 ' and (ds)^ directly 
computed would produce 9j symbols. In fact, so do (us)^ and {ds)(°\ but they produce a 9j 
symbol with one zero argument, which can be written as the product of two 6j symbols, as 



seen below. Analytic forms for 6j symbols with one argument < 2 appear in Edmonds [19 



which is precisely what is needed to compute the matrix elements of tensors up to rank 2. 
The analytic forms of the matrix elements of interest are: 

M (0) = Ss'sSs' z sMsJs'sMs s (-l) 1+S+S ^- Sd 



x ^S U (S U + 1)(2S U + 1)S S (S S + 1)(2S S + 1)(2/' + 1)(2J + 1) 
Sd S u I i I S S x I 



1 r s u 



1 v s K 



(3.18) 



(S Z U ) 



Ss'zSjs^Suds'jSjds'sSs (-1 



^S-S'+Sz+Ss+I'-I-Su-Sj 



X\/S U (S U + 1)(2S U + 1)(2I' + 1)(2J + 1)(2S" + 1)(2S + 1) 

S d S u I )( S s I s)( 1 S' S 
1 V S u 1 S' V f I S z -S z 



(3.19) 



{si) 



SsisJi>i5s>sJs'sMs s (-±) 1+Sz+Ss+ WS s (S s + 1)(2S S + 1)(2S> + 1)(2S + 1) 



I S s S 
1 S' 5, 




(3.20) 



) (2) = ^ 5z ^,^ su ^ sd ^ 5s (-l) 5 - 5 ' + ^ + ^^-^ s - 




(3.21) 



\UU i 

x A /(25 tt -l)5' tt (5 11 + l)(25 11 + 3)(2/' + l)(2/ - \)(2S'+ DC2S + ]) 
S d S u I } \ S s I s)( 2 S' S 

2 s' r\[o S z -S z 

(ss)® = 5 s , s Jri5 SLS Js> d sMss(-l) S ° +Ss+I 

x^f(2S s - l)S a (S a + 1)(2S S + 1)(2S S + 3)(2S' + 1)(2S + 1) 

f I S s S 1 1 2 5" S \ 
x . (3.22) 

[2 S' S s \\0 S z -S z ) 

Note that, in the interest of exhibiting maximal symmetry, the remaining CG coefficients 
have been written as 3j symbols |19 . 



Charge radius and quadrupole transitions are diagonal in both charge and strangeness, 
and thus connect only states of the same values of -Z3 and N s , but do not necessarily conserve 
isospin. Furthermore, charge radius operators are scalars and thus connect only states of 
the same total spin J, but quadrupole operators are rank 2 and therefore can connect spin 
3/2 to 3/2 or 1/2, but not 1/2 to 1/2. These selection rules are reflected by the transition 
matrix elements represented in the tables. Values at arbitrary N c for the matrix elements 
(aP)^ (a 7^ (3) for all relevant states are collected in Table |, for in Table [TI], for (a«)' 2 ' 
in Table |TTT| , and (a/?)*- 2 -* (a 7^ (3) in Table 0. One may obtain results for charge radii 
and quadrupole moments, including arbitrary SU(3) symmetry breaking, by combining the 
results of these tables using the expressions (p.5|)-( pl9|) derived in the previous Subsection. 



IV. RESULTS 



We exhibit in the Appendices expressions for the baryon charge radii (A) and quadrupole 
moments (B) under three sets of assumptions familiar to researchers in the quark model. The 
first is simply to assume no SU(3) flavor symmetry breaking except for that inherent from 
the quark charges. The second and third, which we call "quadratic" (n = 2) and "cubic" 
(n = 3) SU(3) symmetry breaking, respectively, correspond to modifying the spin-spin terms 
in the following way: 

(TiCT j — > (Ti<T jVPb 1 jwbl~ x m^ (4-1) 
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where rrii denotes the constituent mass of quark i, and m = m u = is the light quark 
constituent mass, thus fixing all parameters a-f and b'-f defined in ( |3.2| ) and ( |3.3| ). Let 
us henceforth abbreviate with m/m s = r the degree of SU(3) flavor symmetry breaking. 
The quadratic mass dependence arises in a constituent quark model with the dominant 
interaction mechanism being one-gluon exchange, while the extra mass factor for cubic mass 
dependence arises from a quark propagator between photon absorption and gluon emission. 

One can also mix these pictures so that, for example, the operator labeled by B(B') uses 
cubic SU(3) symmetry breaking, while that labeled by C(C) uses quadratic breaking. Note 
that no SU(3) symmetry breaking has been introduced in the one-body (A) operator. 

We hasten to add that it is not necessary to adopt any model dependence beyond the 
single-photon exchange ansatz, since we have computed and tabulated all the relevant primi- 
tive matrix elements. However, the SU(3)-symmetric, quadratic, and cubic models provide a 
useful picture in which to investigate the consequences of the arbitrary N c calculation, with- 
out adding any new parameters. Another interesting scheme for SU(3) symmetry breaking, 
but not investigated here, is that provided by chiral perturbation theory. In this case, 
loop graphs that include the octet of light mesons as Goldstone bosons produce terms with 
different r dependences, e.g., logarithmic. 

We have also used expressions |2D| for the quark charges that simultaneously guarantee 
chiral anomaly cancellation of the iV c -extended standard model and fix the total charges of 
the iV c -quark baryons to equal their N c = 3 values. These are: 

Qu,c,t = (N c + 1)/2N C , Q dtS , b = {-N c + l)/2N c . (4.2) 

Given the primitive matrix elements of Tables |-fy|, one may alternately compute expressions 
using the strict N c = 3 values for the quark charges, with the caveats that the anomaly 
cancellation conditions are no longer satisfied and the baryon charges become N c dependent. 

One may choose any of a number of schemes for obtaining interesting predictions from 
our results. Rather than selecting just one and producing exhaustive results, we discuss 
several possibilities and exhibit illustrative examples. 

First, one must differentiate predictions that hold in the physical case N c = 3 from those 
that hold in the 1/N C —>■ limit. The former have the advantage of not depending on 
the quality of the 1/N C expansion at N c = 3, but may depend upon delicate cancellations 
between terms at different powers in 1/N C . Thus, such predictions may hold well only for 
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N c = 3 but not for 5,7,. . . . Conversely, the latter have the advantage of holding to a 
desired level of accuracy for small values of 1/N C , but these corrections may turn out to be 
numerically significant for N c = 3. As we exhibit below, a number of relations are found to 
hold independent of N c and thus satisfy both criteria. Moreover, there are even relations 
that hold both for N c = 3 and 1/N C — > 0, but differ in between; we discuss one such example 
with the quadrupole moments below. 



A. Charge Radii 

Baryon mean-square charge radii, denoted here by r 2 B , are defined by the l.h.s. of Eq. ( |3.2| ). 
divided by the total charge if it is nonzero. In all suggested substitutions for the coefficients 
a at ■ ■ ■ i fa, isospin violation is introduced via the single-photon exchange ansatz through a 
single power of the quark charge operator Q a , which transforms as a combination of / = 
and 1=1. Consequently, any combination of charge radii only sensitive to / = 2 or / = 3 
operators must vanish. These are 

1 = 2 : 

= 2r 2 A++ - 3r 2 A+ + 3r 2 A o + r\-, (4.3) 

= r| + -2r| -r|_, (4.4) 

= r|, + -2r|, -r|,_, (4.5) 
1 = 3 : 

= 2r| ++ -r| + -r| -r|_. (4.6) 

Equations Q4.3|) and ( f4.6| ) were first derived in Ref . || , while Eqs. ( fO| ) and ( f4.5| ) are £ equal- 
spacing rules, adjusted for the negative charge of E^~. In addition, using Appendix A one 
finds 3 linear combinations of charge radii with A c -independent coefficients that vanish 
for arbitrary values of N c and r, either in the quadratic or cubic case of SU(3) symmetry 
breaking: 

= -20(r;; + r 2 n ) + 5(r|+ - r|„ + 3r|_) + 5r 2 A - 4(4r^ ++ - r 2 A+ - Ar\ + 7r 2 A -) 

+10(r|. + -r|. + 3r|._), (4.7) 

= 4(r 2 p - hr 2 n ) - (5r| + + 3r| - r|_) + 35ri + 4(2r^ ++ + r 2 A+ + r 2 A0 - r 2 A .) 

-2(5r|.+ +3r|» -r|,_), (4.8) 
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= -8(r; + r 2 ) + (rl+-9r^-) + 2r 2 A + 16r|_ + 2(r^+ - 9r^-) + 32r|,_ - 16r^_.(4.9) 

Since so many of these charge radii are currently unmeasured, these relations are presented 
merely for completeness. 

A number of charge radii can be related if one permits ^-dependent coefficients. A 
particularly pretty example is 

r * = ids* < 4 - 10 ' 

The N c = 3 version of this relation, but applied to magnetic moments rather than charge 
radii, is known from the early days of SU(3) flavor (see, e.g., pi]]). Using the measured 
value for the neutron charge radius fl2.3|) and N c = 3, one predicts r\ = —0.057(3) fm 2 . 



The expression ( 4.10 ) is especially interesting because it definitively predicts — independent 
of SU(3) symmetry breaking — that the sign of the charge radius of A is the same as that 
of the neutron, i.e., negative. This is not the case, for example, in the calculations of 
Ref. |22| (extrapolation from lattice results) or Ref. |23| (constituent quark model). For 



other comparisons, Refs. |23], |24], [25[ also contain tabulated results of various authors' baryon 
charge radii calculations. 

Another interesting result is that 

r| + =r 2 + 0(l-r), (4.11) 

so that the two charge radii are equal in the SU(3) limit for all values of N c , as is evident 
from Eqs. ([All ) and ( A13 ). One expects that this relationship holds to better accuracy than 



r| = r 2 , which requires both the SU(3) limit and N c = 3 to hold, so that both corrections 
of 0(1 - r) and 0(1 - 3/N c ) occur. 

There is a relation of the latter type between the 3 measured charge radii: 

4 ~ rl- + r 2 n = 0(1 - r) + 0(1 - 3/N c ). (4.12) 

Or, if one allows ^-dependent coefficients, 

\{N C - l)(rj - r|_) + rl = 0[(1 - r)(l - 3/i\Q]. (4.13) 

Numerically, the l.h.s. nearly vanishes: One finds —0.06(15) fm 2 . However, the uncertainty 
is the figure of merit here, since it may be used to gauge the typical size of a first-order 
SU(3)-breaking effect. 
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Yet one more interesting result is: 



I- = \ (4 + rl- + r 2 n ) + 0[(1 - rf] + 0(1 - 3/N c ), 



(4.14) 



in both the quadratic and cubic forms of SU(3) symmetry breaking. The second of the order 
terms implies that the result holds exactly only for N c = 3, but this can be pushed to higher 
order if one allows ^-dependent coefficients on the l.h.s., just as in Eq. ( gl3D . The virtue 
of this expression as it stands is that it predicts r§_ up to second order in SU(3) symmetry 
breaking entirely in terms of the known charge radii, with simple coefficients. One finds 



where the uncertainty is dominated by that of the r|_ measurement, which is much larger 
than the piece of the uncertainty one would estimate from the second-order SU(3) symmetry 
breaking. 

Lastly, our expressions for charge radii depend upon only four parameters, A, B, C, and r 
[and a choice of scheme for SU(3) breaking]. Once one additional charge radius is measured, 
it will be posssible to solve for all the parameters and predict all of the other charge radii. 
In this sense, r|_ was chosen in the previous paragraph because it could be predicted in 
terms of just the measured charge radii, with little sensitivity to the SU(3)-breaking scheme 
or the parameter r. Conversely, if it becomes the fourth measured baryon charge radius, it 
will not effectively constrain the parameter r. 

On the other hand, it turns out that the value of (to give one example) cannot be 
predicted in terms of rt n 2 _, even in the SU(3) limit. Turning this unfortunate observation 
around means that a measurement of would provide an extremely sensitive probe of 
SU(3) symmetry breaking, and permit high-quality predictions of all other charge radii. 

Using the three measured charge radii and the assumption of either quadratic (n = 1) 
or cubic (n = 2) SU(3) symmetry breaking and an assumed value of r, one can obtain at 
least ranges for the values of the parameters A, B, and C. Since these parameters have an 
expansion in 1/iV™ starting with m = 0, the values obtained for N c = 3 hold for arbitrary 
N c . One finds: 



r |_ = 0.64(8) fm 2 



(4.15) 



(1 -2r-2r n )(r 2 + r 2 ) + 3r|- 
r l ~ r l- + - r - r n )r 2 n 
r 2 p -rl- -i(l-2r-2r>* 



+2(2 -r 




r n )A, 



r n )B, 



r n )C. 



(4.16) 
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Note that these all become Eq. ( |4.12|) in the SU(3) limit. The uncertainties in these expres- 
sions are dominated by that of r|_, even more than by that of n. Including only the former 
and assuming the value r = 330 MeV/ 540 MeV suggested by the constituent quark model, 
one computes: 

A = +0.56 ±0.29 fm 2 , B = -0.05 ± 0.29 fm 2 , C = -0.33 ± 0.44 fm 2 . (4.17) 

With these central values one predicts r\+ < r 2 for N c = 3, contrary to the physical picture 
that the A is an excited state of the nucleon and hence is more extended in space. However, 
it is well within the error bars of r|_ for B to be substantial and for C to nearly vanish. 
For example, with r|_=0.73 fm 2 as suggested by the upper value of the range of statistical 
uncertainty in Eq. ( |2.4j ) and r as above one obtains A = 0.79, B = 0.18, and C = 0.02, a 
hierachy of parameters leading to the reasonable conclusion r A+ > r 2 



B. Quadrupole Moments 

The operator ( |3.3|) defining quadrupole moments shares with the charge radius operator 
( p.2|) the feature in the single-photon exchange ansatz that isospin violation enters only 
through a single power of the quark charge operator Q a . Thus, as before, combinations 
sensitive only to I = 2 or I = 3 operators must vanish. These are: 



7 = 3 : 

= Q A ++ - 3Q A + + 3Q A o - Qa-, (4.18) 
1 = 2 : 

= Q A ++ - Qa+ - Qa» + Qa-, (4.19) 

= Qa+ p ~ QA°n, (4.20) 

= Q s * + -2Q E ,o+Q s ,-, (4.21) 

= Q s * +S+ - 2Q E ,o s o + Q s .- E _. (4.22) 



The first three of these expressions were obtained in Ref. 0, while the last two are E(*) 
equal-spacing rules, obtained for iV c = 3 in and [^|. In addition, there is precisely one 
linear relation with iV c -independent coefficients that holds for all values of N c in all cases of 
SU(3) symmetry breaking studied here: 

= Q H *- - Q a - - V2Q S *- S - . (4.23) 
15 



Unlike the charge radius case, only the N — > A quadrupole transition matrix element 
has been measured (via photoproduction experiments [0]), and even here the extraction of 
Qn^a is plagued by a large model dependence: 

Q N ^ A = -0.108 ± 0.009 (stat + syst) ± 0.034 (model) fm 2 . (4.24) 

Since there are 3 undetermined parameters (B', C, and r, as well as a choice of SU(3)- 
breaking scheme), we do not attempt to predict any of the other quadrupole moments 
numerically. However, we can still make a number of interesting observations based on the 
structure of the expressions in Appendix B. 

First note that each quadrupole moment expression is 0(N®) or smaller in the 1/N C 
expansion. Indeed, only one coefficient, B', contributes to this leading order (in the case 
of charge radii, both A [for charged baryons only] and B contribute at leading order). A 
few moments' study will confirm that the diagonal quadrupole moments are given by the 
expression 

Q(I 3 , Y) = h[l - y(4 - Y)/lh\B' + 0(1/N C ). (4.25) 

Proportionality to 73 also holds for the leading terms of the transition quadrupole moments 
when the initial and final baryon states have the same value of I (£*£ and 5*H). This 
behavior, due to the dominance of the isovector portion of the quadrupole operator, is 
familiar from Skyrme model and a variety of other model calculations |2IJ . Taken at face 
value, it predicts an appreciable quadrupole moment for A , as well as Qq- = 0. However, 
the subleading corrections in 1/N C , which interpolate between the extreme N c — > oo and 
N c = 3 cases, soften this behavior. In particular, the diagonal quadrupole moments in the 
strict N c = 3 case with no SU(3) flavor symmetry breaking obey 

Q = 4g/3 (B' + C/3) , (4.26) 

where q is the baryon charge. 

As discussed in Sec. |T|, the charge radius and quadrupole operators are very similar in 
that both represent spin-dependent electromagnetic couplings to baryons. It thus should 
not be surprising that their coefficients are related in a typical model. In Ref. we saw 
that the one-gluon exchange picture gives rise to the relation 



1 2 N c N c + 5 

Q ^ = ^J^T-^—r (427) 
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for which the factor on the r.h.s. following r\ equals unity both for N c = 3 and N c — > oo. In 
fact, this result can be obtained in a much more general setting. The key constraints needed 
to obtain Eq. ([£270 are B' = -2B, C = -2C, and as argued in Ref. JT§| for N c = 3, the 
same relation may be derived using not only one-gluon exchange, but one-pion exchange or 
scalar exchanges, or a mixture of these. 

A particularly useful measurement for determining the coefficients in the quadrupole 
sector would be that of Qn-- As we see in Appendix B, the value of Qn- is very sensitive 
to the precise nature of SU(3) symmetry breaking, even more so than r^_. We discussed 
in Ref. H ideas in the literature for the experimental determination of Qn-', while such 
experiments are challenging, they appear to be feasible. The value of Qn- would teach us 
much about the shape of baryon wave functions and the nature of SU(3) flavor symmetry. 



V. CONCLUSIONS 



We have presented techniques that permit the calculation of matrix elements of operators 
acting upon 3-flavor baryon states with arbtitrary N c . The approach requires only SU(2) 
Clebsch-Gordan coefficients and combinations of them in the form of 6j symbols. We tabu- 
lated the values of a set of primitive operators for all relevant states and demonstrated how 
they can be combined to give results for charge radii or quadrupole moments, using any 
chosen pattern of SU(3) (or isospin) flavor symmetry breaking. In particular, we presented 
in the Appendices results using the single-photon exchange ansatz, augmented by either no 
SU(3) symmetry breaking, or one of two popular types of SU(3) symmetry breaking sug- 
gested by (but not limited to) the quark model with gluon exchange. We obtained a large 
number of interesting predictions and demonstrated how many others can be made either 
working in some particular model, or once a small number of additional baryon charge radii 
or quadrupole moments are experimentally measured. 
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APPENDIX A: CHARGE RADIUS EXPRESSIONS 

We use the subscripts 0, Q, and C to denote expressions with zero, quadratic, and cubic 
SU(3) flavor symmetry breaking via constituent quark masses [see Eq. (|4.1|) 1 in addition 
to that provided by the quark charge operator. In the SU(3) symmetry limit (r = 1) the 
expressions with subscripts Q and C reduce to those with subscript 0. 
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APPENDIX B: QUADRUPOLE MOMENT EXPRESSIONS 



We use the subscripts 0, Q, and C to denote expressions with zero, quadratic, and cubic 
SU(3) breaking via constituent quark masses [see Eq. (|4.1]) 1 in addition to that provided 
by the quark charge operator. In the SU(3) symmetry limit (r = 1) the expressions with 
subscripts Q and C reduce to those with subscript 0. 
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TABLE I: Values of N uAs [whence (S 2 a ) = (aa)^ = (A Q /2)(A a /2 + 1)] and matrix elements of 
the rank-0 tensors {af3)^ with a / 0. Since spin is unchanged by these operators, the matrix 
elements vanish for all off-diagonal transitions except S°A; in that case, the only nonvanishing 
entries are («s)(°) = -(ds)^ = -§ y/(N c - 1)(A C + 3). 



State 




Aw 


A 






(ds)W 
\ / 


A++ 
ZA 


1 f TV i "J 1 ! 

2 V-/V c + 6) 


1 / 7\r q\ 
2(iV c - 6) 


n 
u 


1 f 7V AT _l_ 7"l 


n 
u 


n 
u 


A + 


1 f TV i 1^ 


1 ( AT 1} 


n 
u 


1 T" 7\/"2 I /I AT" OQ\ 


n 
u 


n 
u 


AO 




1 / 7\r i i ^ 


n 
U 


1 1 7\r2 l /I AT" OQ\ 


n 
U 


U 


A — 
ZA 


1 1 AT 


1 ( AT i 


n 
U 


1 1 AT TV _l_ 7\ 


n 
U 


n 
U 


V* + 
Zj 


1 f TV ±11 


1 / AT" q\ 


1 


1 1 AT AT _1_ K\ 


i i Tat _i_ ^\ 


1 z' AT ^ \ 


v*o 

Zj 


1 ( AT -\ \ 


1 ( AT -\ \ 


1 


1 ( Aj2 i HI 1 n\ 


,1 
+ 4 


,1 
+ 4 


v^* — 

Zj 


1 / AT 0\ 


1 / 7\7" _i_ -1 \ 


1 




1 ( AT "i\ 


i 1 / AT" _j_ r \ 




1 ( AT 1} 


1 1 AT OA 


Z 


1 ( AT Q\ / 7\7" i q\ 


I 1 / AT i o\ 
+ 12 V iV c + <3j 


1 / AT" q^i 




1 / 7V7" q\ 


1 / A/" 1 \ 

2l-< v c — 1J 


o 
Z 


1 1 AT QU A7" I 9\ 


1 f' AT 1} 


1 1 / AT _J_ Q\ 

+ 12 V iV c + o) 


o 

1 i 


1 / AT <1\ 

2{I\ C - 6) 


1 / AT <1\ 


q 


1 / AT" Q\ / AT" I 1 


n 
u 


n 
u 


p 


1 ( AT i 1 \ 

2-(iV c + 1) 


1 1 AT 1 

2(iV c - 1) 


u 


1 / AT 1 \ / A7" i cr \ 

-Je( N c - l)(A c + 5) 


U 


U 


n 


2(JVc-l) 


i(A c + l) 





-^(A c -l)(A c + 5) 








£+ 


^(iVc + 1) 


^(A c -3) 


1 


-^(A c -3)(A c + 5) 


4(^ + 5) 


+1(^0-3) 


S° 






1 


-^(A C 2 + 2A C -19) 


l 

2 


l 

2 


A 


5(JVc-l) 


5(^c-l) 


1 


-^(A c -l)(A c + 3) 








£- 


i(A c -3) 


i(A c + l) 


1 


-^(A c -3)(A c + 5) 


+1(^-3) 


-1(^0 + 5) 


S° 


5(JVc-l) 


i(A c -3) 


2 


-B(^e-3)(A c + 3) 




+H(^c-3) 




|(JVc-3) 


5(^Vc-l) 


2 


-^(A c -3)(A c + 3) 




-|(^c + 3) 



26 



TABLE II: Matrix elements of the operators S*, 5J, and S% in the state of maximal S z . 
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TABLE III: The matrix elements (aa)' 2 ' for all relevant states. 
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TABLE IV: The matrix elements (a(3)( 2 \ with a//?, for all relevant states. 
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